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Abstract 

We modify the Schrodinger equation in a way that preserves its main properties but makes 
use of higher order derivative terms. Although the modification represents an analogy to the 
Doebner-Goldin modification, it can differ from it quite distinctively. A particular model of 
this modification including derivatives up to the fourth order is examined in greater detail. We 
observe that a special variant of this model partially retains the linear superposition principle 
for the wave packets of standard quantum mechanics remain solutions to it. It is a peculiarity 
of this variant that a periodic structure emerges naturally from its equations. As a result, a free 
particle, in addition to a plane wave solution, can possess band solutions. It is argued that this 
can give rise to well-focused particle trajectories. Owing to this peculiarity, when interpreted 
outside quantum theory, the equations of this modification could also be used to model pattern 
formation phenomena. 
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Recently there has been a considerable interest in the Doebner-Goldin (DG) modification of the 
Schrodinger equation [|TJ. (See also |2j for a more complete list of references to this subject and 0, f|] 
for the latest update on the progress in a broader context related to the matter in question.) It is the 
purpose of the present letter to propose yet another modification of this fundamental equation that, 
similarly as the Doebner-Goldin modification, makes use of the current formulation and preserves the 
main features of the equation discussed, these of homogeneity || and weak separability of composed 
systems 0. Moreover, a subset of equations of this modification, defined by certain values of its free 
parameters, complies with the Galilean invariance, giving rise to a special version of it that, except 
for linearity, possesses all the main properties of the Schrodinger equation deemed physically relevant. 

Let us note that even though in some cases the homogeneity of nonlinear generalizations of the 
Schrodinger equation does entails their weak separability |7|], the properties mentioned are, in general, 
independent ||. As pointed out in ||, the homogeneity of nonlinear variants of this equation is essen- 
tial for a unique defintion of their energy functionals. This property is also necessary for modifications 



of the Schrodinger equation to possess a weakly separable multi-particle extension [g, |T(J . 

What differs our proposal from the DG modification is the use of terms that involve derivatives 
of the order higher than second and of higher polynomial degrees. As we will see, the simplest 
extension of this kind employes derivatives of fourth order and of second degree. The higher degrees 
are required for the completness of the formulation. As we will also see, a certain, physically most 
attractive, variant of this modification should be viewed as an extension of the Schrodinger equation 
rather than its modification for it changes the main properties of this equation less dramatically than 
the majority of the modifications. In particular, this true about its solutions. 

To begin with, let us first reformulate the DG modification in a way that will be convenient for 
the intended generalization. To this end, let us write the modified Schrodinger equation as 

where 

n 

F {x} =5>iFi 

i=l 

and Xi are some dimensionless coefficients that form a generic array {x} while [\&, \&*] are functionals 
of \I/ and homogeneous of degree zero in these functions. The coupling constant D has the 
dimensions of the diffusion coefficient, meter 2 second _1 , in the DG modification. In what follows, 
we will work with p = R 2 = \ff\ff* and S, the probability density and the phase of the wave function 
^ = Rexp(iS), correspondingly. The general form of the functional employed by Doebner and Goldin 
is 

F ^f[p j S]=x 1 AS + x 2 VS-(^\+x 3 ^+x A (^A +x 5 (ys)\ (2) 

The imaginary part of the Schrodinger equation leads to the continuity equation. The standard way 
to obtain it is to multiply both sides of the Schrodinger equation by \&* and take the imaginary part 
of the ensuing expression. The result turns out to be 

dp + Ay • (pVS) = (3) 



dt m 

with the probability current identified as j = hpVS/m. 
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In the same manner, we would like pF^ a y to form the divergence of some current. One can show 
that two terms emerge to play this role: V • (pVS 1 ) and Ap. One obtains these in a unique way by 
putting di = a<i = a and a 4 = = 0. Renaming a^D — ► D' and a/} — ► D allows us to write the 
modified continuity equation asQ 

% + -V • (pVS) + ■ (pVS) + D'Ap = (4) 

or m v 7 v 7 

which can also be put in the form 

t + ^V-(pV5) + /rAp = 0, (5) 

where m* = mj (1 + Dm/h) = m//3 is the effective mass of a quantum system. Let us notice that the 
most uniform way to write the last equation is as follows: 

| + Av. (p v S) + C 'V.(^)=0. ,6, 

If we further note that S = | In ( 1 i r */ 1 i r ) and Vp/ p = In 1 i r * 1 i r , we see that this way reveals the ubiquitous 
role of the logarithmic function in the DG modification^] and suggests that the probability current 
be considered as consisting of two components having velocities vg = ^VS* and v p = D'Vp/ p. The 
latter is the diffusion component. 

The linear Schrodinger equation is invariant under the Galilean transformation of coordinates 
x = x + vt,t = t provided the phase of the wave function transforms as 

S(x, t) — S'(x , t') + mv ■ x + -mv 2 t . (7) 

One observes that the DV ■ ^pV5*^ component of the additional current would break the Galilean 
invariance of the continuity equation if it were not for the mass redefinition. To ensure the Galilean 
invariance of this equation it is also required that the effective mass replaces the "bare" mass every- 
where a reference to the latter is made, including in particular the phase of a quantum system in the 
above transformational formula. However, the redefinition in question by no means guarantees the 
Galilean invariance of the entire Schrodinger equation, which is maintained only if D = 0. 

It is also easy to notice that a special version of DG modification can be straightforwardly derived 
by utilizing VS 1 and Vp/p as the components of a general vector potential Adg = d\VS + g^V ' p/ p. 
In the presence of Adq the free Schrodinger equation reads 

One arrives from it at the linearizable variant of the DG proposal that generates the continuity equation 
in the form (4) or (5) and modifies the other part of the Schrodinger equation by adding terms (VS") , 

1 In the original DG formulation |Q, D does not make its appearance at all. The case of nonzero D is treated in a 
subsequent generalization of their scheme involving nonlinear gauge transformations || . 

2 This should not come as a surprise as the logarithmic function was already used for the sole purpose of ensuring 
that the modification of the Schrodinger equation proposed in B obeys the weak separability condition. 
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(Vp/p^j , and VS ■ Vp/p. The coefficients that stand by these terms depend only on two constants 
which are functions of di. It is clear now why this variant is linearizable and how to transform it 
into the linear Schrodinger equation. The latter is achieved by a nonlinear gauge transformation ||, 
similarly as this is done in electrodynamics. Some caution should be exercised when dealing with the 
VS part of the gauge potential or the corresponding nonlinear gauge transformation. As pointed out 
above this part is responsible for the mass redefinition in the continuity equation. However, instead 
of the mass redefinition, one can also enforce the Galilean invariance of this equation by redefining 
the phase S' = (3S. Now, we see that this alters the range of the phase from 2ir to 2irf3 and should 
be taken into account when dealing with the uniqueness of the wave function Yet another way 
to ensure the Galilean invariance is to redefine the Planck constant, h — > f3h, which eventually brings 
the nonlinear Schrodinger equation to the form of its linear prototype, except that with the rescaled 
Planck constant. This completes the linearization process. Being linearizable, the discussed variant 
of the DG modification is physically equivalent to the linear Schrodinger equation, and, in particular, 
it satisfies the Ehrenfest equations. The DG proposal in its full form is not linearizable. However, 
similarly as the variant in question is related to the linear Schrodinger equation via nonlinear gauge 
transformations and thus constitutes a kind of family with it, the equations of the full DG modification 
form a family that is closed under more general gauge transformations 0. 

No stipulation is made on the coefficients {b}. However, if one wants the modified Schrodinger 
equation to be Galilean invariant, (7) implies the unique choice of 62 = &5 = 0. A restricted version 
of the DG modification for b 2 = and 26 3 — 64 = in F^y ty*] was shown || to be derivable from 
a local Lagrangian density. 

The main idea of our proposal is to extend the DG modification so as to include higher derivative 
terms. Despite an apparent simplicity of this task, the resulting construction turns out to be rather 
complex and can have diametrically different features than the proposal of Doebner and Goldin. 
We will first consider the leading order case, involving the derivatives up to the fourth order and 
the second degree, which, in addition to its simplicity, seems to be the most promising on physical 
grounds, and then suggest how to generalize this construction to allow also higher order terms. It is 
rather straightforward to convince oneself that 



rpext 
*{*} 




[p,S] = x 1 AAS + x 2 a(^J + x 3 A +X4 a^-VsJ + 

X5 A (vsY + ■ V (AS) +X7 fl.v(^) + 

' p p V p J 

X9 ^.v(^.vs] +Xl0 ^.v(vs) 2 + 
p \ p J p v ' 

x n VS ■ VAS + x 12 VS ■ V {^j j +xi 3 VS -V {^j ( 9 ) 

is the right choice for the homogeneous terms up to the fourth order. Terms (AS) 2 , AS (Ap/p), and 

AS (Vp/p^j have not been included. Even if homogeneous of degree zero in both \I/ and \^*, they 
would nevertheless break the weak separability of the modified Schrodinger equation. 

3 See also |fTTf for a more profound discussion of this issue in the case the nonlinear gauge transformations are 
supposed to form a group. 
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To cast more light on this issue, let us demonstrate the weak separability of the Schrodinger 
equation in the hydrodynamic formulation. We are considering a quantum system made up of two 
noninteracting subsystems in the sense that || 

V(x 1 ,x 2 ,t) = V 1 (x 1 ,t) + V 2 (x 2 ,t). (10) 

We will show that a solution of the Schrodinger equation for this system can be put in the form 
of the product of wave functions for individual subsystems for any t > 0, that is, *ff(xi,x 2 ,t) = 
t)^/ 2 (x 2 , t) = Ri(xi, t)R 2 (x 2 , t)exp {i(Si(xi,t) + S 2 (x 2 , t))} and that this form entails the sepa- 
rability of the subsystems. The essential element here is that the subsystems are initially uncorrelated 
which is expressed by the fact that the total wave function is the product of *&i(xi, t) and ^ 2 (x 2 , t) at 
t = 0. What we will show then is that the subsystems remain uncorrelated during the evolution and 
that, at the same time, they also remain separated. It is the additive form of the total potential that 
guarantees that no interaction between the subsystems occurs, ensuring that they remain uncorrelated 
during the evolution. However, such an interaction may, in principle, occur in nonlinear modifications 
of the Schrodinger equation even if the form of the potential itself does not imply that. This is due to 
a coupling that a nonlinear term usually causes between ^i(xi,t) and ^ 2 (x 2 ,t). As a result, even in 
the absence of forces the very existence of one of the particles affects the evolution of the other one, 
clearly violating causality. 

The discussed separability is called the weak separability since it assumes that the wave function 
of the total system is the product of the wave functions of its subsystems in contradistinction to 
the strong version of separability that does not employ this assumption. As shown by Liicke [O, ITS 



weakly separable modifications, such as the modification of Bialynicki-Birula || or the Doebner-Goldin 
modification, |]] can still violate separability when the compound wave function is not factorizable, and 
thus they are not strongly separable. An alternative effective approach to the strong separability has 



been proposed by Czachor JL4]. This approach treats the density matrix as the basic object subjected 
to the quantum equations of motion which are modified^ compared to a nonlinear Schrodinger equation 
for the pure state. It admits a large class of nonlinear modifications including those ruled out by the 
fundamentalist approach advocated by Liicke and even those that are not weakly separable as, for 
instance, the cubic nonlinear Schrodinger equation. 

The Schrodinger equation for the total system, assuming that the subsystems have the same mass 
m, reads now 



dR 2 



,dR 2 



dt 

h 2 h 2 
+—R 2 2 V! ■ (R 2 VxSi) + —R\V 2 ■ (R 2 2 V 2 S 2 ) = R\R 



m 

*± 
m R{ 



dt 1 dt 

' 1 dR 2 



R\ dt 



2 Vi 



(RlVxS, 



+ 



' 1 dR 2 
% 



+ 



h 2 l 



R\ dt mR\ 



V 2 -[RlV 2 S 2 



'IV, 



and 



4 What this means in practice is that the basic equation is the nonlinear von Neumann equation |iq] instead of some 
nonlinear Schrodinger equation for the pure state. 
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— (Ai + A 2 ) RiRz - 2hR l R 2 d ^ Sl ^ S ^ - —R X R 2 (v l S 1 + V 2 S 2 ) 2 - 
m ot m K ' 

h 2 h 2 dS OS 

(Fi + V 2 ) R1R2 = —R2&1R1 + —RiA 2 R 2 - 2hR 1 R 2 —± - 2hR 1 R 2 —l 
mm ot ot 

;2 



+—RxR 2 (Vi^) 2 + — R 1 R 2 (V 2 S 2 Y - VxRiR 2 - V 2 R X R 2 



f 
m 

R\R 2 
h 2 



m 



h 2 A X R X 
m Ri 

(v 2 s 2 ) 2 -v 2 



2h ^- + - 

ot m v 



Vi 



r h 2 A 2 R 2 
m R 2 



+ 0S 2 
2h^+ 



Ot 



'12) 



where we used the fact that Vi ■ V 2 = 0. What we have obtained is a system of two equations, each 
consisting of terms (in square brackets) that pertain to only one of the subsystems. By dividing the 
first equation by R\R 2 and the second one by R\R 2 , one completes the separation of the Schrodinger 
equation for the compound system into the equations for the subsystems. Moreover, we have also 
showed that indeed the product of wave functions of the subsystems evolves as the wave function 
of the total system. A similar analysis applied to (9) convinces us that the chosen functional does 
have the required property of weak separability. Of course, as already mentioned, this by no means 
guarantees that the separability will be maintained for nonfactorizable compound wave functions. 
However, the proposed modification should be strongly separable in the effective approach [pl[| . 
The coefficient D has now the dimensions meter 4 second _1 . By choosing 

and an — a i2 — a i3 — one obtains the following continuity equation 



% + -V ■ ( P VS) + D 1 V ■ (pVAS) + D 2 V 

Ot m v ' v ' 



:«(?) 





















pv (vs) 2 













+ 



(13) 



where D; = a^D. The new currents are revealed to be 



, J4 



Ji = D lP VAS, j 2 = D 2 pV (^j , j 3 = D- iP V 

j 5 = D 5P V(VS) 2 . 
The continuity equation (13) can also be written in the form 

% + ^-( P VS)+DV.[pVF°f[ P ,S] 



zvv ( ^ ■ vs) , 



0. 



(14) 



(15) 



Similarly as in the DG modification, no condition is put on the coefficients {6}, but if one requires the 
modification to be Galilean invariant one should demand that 64 = 65 = bg = 610 = 611 = b\ 2 = 613 = 0. 
For the continuity equation to be Galilean invariant one needs to set -D 4 = D 5 = 0. 
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One can ask if using VFRf [p, S] as a vector potential in the linear Schrodinger equation would 
generate a special variant of the proposed modification similarly as it was demonstrated for the 
linearizable part of the DG equations. It turns out that this would not be so, which indicates that 
the modification discussed may not be linearizable and therefore may contain some new physics that 
cannot be described by the standard quantum theory. Such a formulation would contribute the square 
of the vector potential to the real part of the Schrodinger equation. However, none of its terms, as for 

instance , belongs to the modification under study, being either of the higher order or degree 

in derivatives. Therefore, to accomplish this, an extention of our proposal that would incorporate 
such terms into it would be necessary. Even though this seems to be rather a natural generalization 
of our work, it is not the purpose of the present paper, and because of that it will not be discussed 
here. 

It is not clear if there exists a local Lagrangian for the modification proposed. It should be noted 
that not all equations of interest for mathematical physics are derivable from local Lagrangians, the 
best case in point being the celebrated Navier- Stokes equations with which no such Lagrangian can 



be associated fll6[ . It is also still unknown if the DG modification in its fully developed form can be 
derived from any local Lagrangian density. 

An essential feature that differs this proposal from the DG modification is the existence of un- 
modified stationary states of the Schrodinger equation, characterized by S = —Et/h, for a certain set 
of parameters. For this to occur, one needs to choose D 2 = D 3 = and 62 = b 3 = 67 = 6 8 = 0. If 
in addition to that we assume the Galilean invariance, the Hamiltonian of the nonlinear Schrodinger 
equation reflecting these constraints reads 



iftV' (pVAS) T/ _ 

Hue = H l + H nl = H l + —± ^ '- + hb x AAS + hb 6 -^-V (AS) , (16) 

2 p P 

where Hi denotes the Hamiltonian of the linear Schrodinger equation.^ There may exist other sta- 
tionary states as well. The energy of stationary states is defined as the expectation value of Hme 

Eme = J dVV*H ME y. (17) 

It is straightforward to show that this definition leads to 

Eme = E L + h{b 1 - b s ) J dVpAAS, (18) 

where we assumed that the total current vanishes on a boundary in the infinity and where El is the 
expectation value of H L given by 



+ VR 2 \ . (19) 



The general condition for stationary states, dp/dt = 0, leads to the equation 

— V • (pVS) + D l V • (pVAS) = £>iV • [p (WAS + wVS)] = 0, (20) 



5 The subscript ME stands for the minimal higher order extension of the Schrodinger equation as explained 
further. 
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where u = -ft—, which should be treated along with 



ni 



Ap V2 _ 2 (rSl + V + — (vsY + ^AAs] p 1 ' 2 - hb^P- • V (AS) p 1 ' 2 = 0. (21) 
\ at 2m v ' 2 J p 



Together, they constitute a system of equations from which new stationary states can be obtained. It 
is conceivable that these equations are satisfied by wave functions whose phase S depends also on the 
position. If D\ > 0, one observes that even in the absence of the potential these equations imply the 
existence of a periodic structure with a period 2ir / y/u that emerges naturally as a harmonic solution 
for VS. It is a peculiarity of the discussed model that, unlike in some other physical situations, this 
structure is not assumed but appears as a consequence of the equations of motion. A free particle, 
in addition to a plane wave solution, would display solutions characteristic of a particle moving in a 
crystal, the band solutions. As observed in [ 17] , an electron in a crystal is well focused in the same 



sense a particle in an accelerator is. The equation that governs strong focusing in accelerators is the 
same as the one that describes the motion of the electron in the periodic crystalline lattice and entails 
the band structure of solids. It is the Mathieu equation. In the discussed situation, one obtains 
this equation from (21) by putting b§ = 0. The strong focusing that this equation implies may offer 
a novel solution to the problem of particle trajectory in quantum mechanics. As is well known, in 
standard quantum mechanics the particles are represented by wave packets. Yet, the wave packets in 
linear theory spread without limit, and in so doing they contradict the existence of sharp well defined 
particle trajectories that one observes in, say, bubble chambers. In nonlinear quantum mechanics this 
problem can, in principle, be solved. This is made possible by particle-like solutions, i.e., localized 
and non-spreading configurations known as solitons that are generic for nonlinear equations of motion. 
Being well localized, the solitons produce sharp trajectories in the bubble chamber. But, as we noted, 
the equations of the modification under study can possibly give rise to sharp well focused trajectories 
even if the particles are not solitons. Interestingly, it is the phase of the wave function that entails 
this remarkable solution, adding a new twist to the meaning of the wave-particle duality. 

On the other hand, if treated outside the quantum theory, the discussed periodic solutions might 
provide a viable model for the pattern formation phenomena |18[. It should be noted in connection 



with this that the equations of the Doebner-Goldin modification have their predecessor in a phe- 
nomenological equation which originates in the study of the wave propagation in fluids and plasmas 



with sharp boundaries and dissipation [19 



As we see, with the assumptions of Galilean invariance and existence of unmodified stationary 
states, the number of undetermined constants of our modification reduces to only three, D\ and bi, 
and be- This particular variant of the modification supports ordinary Gaussian wave packets for which 

c mtx2 1 4- to foo\ 
S = — tt, nT arctan — , (22) 

2 (t 2 + t§) 2 t ' V ; 

(for simplicity in one dimension and in natural units) and the coherent states satisfying AS = 0, and 
this again distinguishes it not only from the DG modification that does not allow for any of these 
packets, but also from the modifications of Staruszkiewicz |2l] and the one proposed in |§ which 



admit only the coherent states, excluding the ordinary wave packets. The fact that the discussed 
variant admits ordinary wave packets indicates that its nonlinearity is weak. It is a general property 
of nonlinear modifications of the Schrodinger equation to exclude such packets. Since both ordinary 
Gaussian wave packets and coherent states constitute the result of superposing more elementary wave 
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functions, each of which is a solution to the modification concerned, this means that the variant in 
question maintains, if only partially, the linear superposition principle. We choose to call this novel and 
rare property among nonlinear modifications of the Schrodinger equation the weak nonlinearity. To the 
best of our knowledge, the property in question is shared by only two other nonlinear modification of 
this equation [[H], ^0[. We will call this variant of the modification the minimal higher order extension 



of the Schrodinger equation for it departs from this equation in the most minimal way, preserving all 
of its standard properties, including the stationary solutions. 

One can easily show that the discussed version of the modification can also be put in the form 
that involves a vector potential A = aVAS as follows 



i— = V-zA * A 2 <Z> + — V ■ + — c 2 Re(A ■ W) + 2Im(A ■ W) . (23) 

dt 2m \ ' 2m 2m 2m L - 1 

The constant a = — m ? x has dimensions meter 2 , the other constants are defined as ci = and 

_ 2mb 6 
° 2 a ■ 

In general, the only notable exception to this rule being the modification of Bialynicki-Birula and 
Mycielski, nonlinear modifications of the Schrodinger equation do not have the classical limit in the 
sense of the Ehrenfest theorem. The discussed variant of the modification is one of such cases. The 
nonlinear terms it introduces entail corrections to the Ehrenfest relations. We will now work out these 
corrections. For a general observable A one finds that 

where the nonlinear contribution is due to H^l = Hr + iHj, Hr and Hi representing the real 
and imaginary part of Hnl, respectively. The brackets <> denote the mean value of the quantity 
embraced. Specifying A for the position and momentum operators, one obtains the general form of 
the modified Ehrenfest relations || 

mj t (r) = (p) + h, (25) 
j t <P> = - (W) + h, (26) 



where 



h = X / dvf P Hl ' ( 2? ) 

I 2 = J dVp (2HiVS - VH R ) . (28) 

In the derivation of the last formula it was assumed that / dW (pHi) = which indeed stems from 
the continuity equation. For the extension in question these integrals are found to be 

I 1ME = D 1 m J dVrV ■ (pVAS) , (29) 



huE = h J dV 



DiV ■ [pVAS] VS - b lP VA 2 S - h P V ( ^ • WAS 



(30) 



One can straightforwardly generalize this construction to even higher order derivatives. To this 
end, we note that the total current functional in (15) is what was the functional from which to build 
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the DG modification. Similarly, the next order current functional will be the functional of formula 
(9). The complexity of the higher order extensions makes them difficult to study. On the other hand, 
this may even not be physically justifiable or interesting. Equations involving fourth order derivatives 
do occur in nonlinear equations aimed at modeling physical phenomena, such as, for instance, the 



pattern formation |]I8 |. Higher order derivatives are rather uncommon. 

The question that deserves further study concerns the new currents. It is quite natural to inquire 
if they form any algebraic structure. Of particular interest is the issue of strong separability of the 
modification proposed in the fundamentalist approach. 

It would also be interesting to compare this modification with the one put forward by Staruszkiewicz 
pl| and extended by this author |L0| which is neither homogeneous nor possesses the property of weak 
separability and with the modification proposed in || which is homogeneous but does not automat- 
ically admit the weak separability of composed systems. The common feature of these modifications 
is the presence of higher degrees of derivatives or derivatives of the order higher than second. Similar 
comparisons would not be possible for modifications with lower degrees of derivatives. This provides 
a unique framework for a better understanding of the concepts of weak separability and homogeneity, 
their physical impact and possible connections. 

The progress in these matters will be reported elsewhere. 
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